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Abstract 

The most general large iV eigenvalues distribution for the one matrix model is shown 
to consist of tree-like structures in the complex plane. For the m = 2 critical point, 
such a solution describes the strong coupling phase of 2d quantum gravity (c = non- 
critical string). It is obtained by taking combinations of complex contours in the matrix 
integral, and the relative weight of the contours is identified with the non-perturbative 
"^-parameter" that fixes uniquely the solution of the string equation (Painleve I). This 
allows to recover by instanton methods results on the non-perturbative effects obtained by 
the Isomonodromic Deformation Method, and to construct for each ^-vacuum the observ- 
ables (the loop correlation functions) which satisfy the loop equations. The breakdown of 
analyticity of the large iV solution is related to the existence of poles for the loop operators. 
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The discovery of the "double scaling solutions" of the matrix models |1J @ led to 
important progress in the understanding of string theories in d < 2 backgrounds and of 2d 
gravity (see [|J for reviews). However, the important issue of the non-perturbative status 
of some of these theories remains unclear, in particular for 2d gravity coupled to unitary 
matter for c < 1. In this letter, we discuss some of these questions in the framework of the 
Hermitian one matrix models. We shall show that a simple generalization of the complex 
integration contour prescription || , which allows to construct non-perturbative — but 
in general complex — solutions of the string equations and of the continuous loop equations, 
leads to real non-perturbative solutions of these equations. This generalization, which 
consists in taking combinations of inequivalent integration contours, has been already 
discussed by Fokas, Its and Kitaev H in the framework of the Isomonodromic Deformation 
Method (IDM) approach to the string equations || , but does not seem to have attracted 
much attention. Our treatment is based on the BIPZ solution of the one matrix model 
ToU , and follows our previous analysis of [11|. We shall show that in the limit N —>■ oo, 



new solutions for the eigenvalues (e.v.) distribution exist, which have not been discussed 
before. They correspond to a distribution of e.v. along tree-like structures in the complex 
plane. Moreover, these solutions depend non-analytically of the coupling constant of the 
matrix model, and will be associated with the sectors with an infinite number of poles 
of the string equation solutions. The non-perturbative parameter which characterizes 
the non-perturbative solutions is simply related to the different weights chosen for the 
contours, and our treatment allows to recover easily by instanton methods some results of 
]7|]||. In addition, we show that to each real solution of the string equation is associated a 
prescription for the asymptotics of the loop operators which defines uniquely observables 
(i.e. macroscopic loop v. e.v.) which obey the loop equations. Finally we shall show that 
these new solutions allow to explain the properties of the solutions for the double well 



matrix models recently discussed by Brower, Deo, Jain and Tan |T2 



In the matrix model formulation of 2d gravity, the partition function F (sum over 
orientable connected 2-dimensional Riemannian spaces) is discretized into a sum over tri- 
angulations, and is written as the logarithm of the partition function Z for the Hermitian 
one-matrix model (F = In Z) , which after diagonalization of the matrix $ can be written 
as an e.v. integral 

r N 

Z N = C N n^^ A -W 2 5 A -( A *) = II( A *- A i) (!) 

J i=l i<j 

where Cat is a normalization factor, the Vandermonde determinant and V the matrix 
potential. The integral ([!]) can be calculated in terms of the matrix elements of the 
operator Q : 7r n (A) i— > Xn n (X), where the 7r n are orthonormal polynomials for the measure 
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dXe NV W . In the double scaling limit, iV — > oo and V — > V^iticai while x = 1 — becomes 
a continuous parameter. Then Q becomes a second order differential operator of the form 

Q = -^2+ M*) (2) 

u is the string susceptibility 

d 2 F 

u(t) = (3) 

where t is the renormalized cosmo logical constant. For the m = 2 critical point (pure c = 
gravity), t scales with iV as t ~ iV 4 / 5 , and u satisfies the Painleve I string equation 

1 d 2 u 



6 dt 2 



+ u 2 = t ; u ~ Vt t -> +oo (4) 



It is known that (|]) fixes uniquely the terms of the asymptotic expansion of u (in powers 
of £( 1-5fc )/ 2 ) as t — > +oo, but that the corresponding series is not Borel summable, and 
that the solutions of (|j) form a one parameter family of "simply truncated solutions" |13|] , 
which differ by exponentially small terms of the form 

8U OC t-l/8 e -|2V3t 5 / 4 (5) 

The real solutions of (|]) have an infinite series of double poles (with residues 1) on the 
negative real axis. If one divides the complex plane into 5 sectors s = I, . . . , V (which 
correspond to (s — 1) 2ir/5 < Arg(t) < s27r/5), these poles extend to an infinite network 
of poles in the sectors H, HI and IV, so that the asymptotics u ~ y/t holds only in the two 
pole-free sectors I and V. 



It was suggested in |5J, and shown more precisely in |TI[]0[§[, that, if one constructs 
the m = 2 theory by starting from a cubic potential of the form 

V(X) = -A 3 + ... (6) 

and defines the integral (|TD by choosing as integration contour for the A^'s the complex 
contour C + (resp. C_) which goes from — oo to joo (j = e l7r / 3 ) (resp. joo), one obtains 
the "simply truncated solution" u + (resp. «_) of (f|) which has poles only in the sector H 
(resp. IV) and satisfies the asymptotics u ~ \fi in the remaining four pole-free sectors . 
From these two solutions, which are analytic on the real axis, one can construct without 
ambiguity the operator Q (0), which is not Hermitian anymore, and the loop operators 



w(p) [0 (where p is the loop momentum), which satisfy the loop equations [p[||lql [16] 



In fact a straightforward generalization of this prescription is to consider linear com- 
binations of the two contours, i.e. to replace 



dX — ► c + / dX + c_ / dX (7) 
c± Jc+ 
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Indeed, the partition function Z will still be real (for real V) if the weight c± are complex 
conjugate 

c± = l±iB (8) 

With this prescription, the orthogonal polynomial method still works, and the recurrence 
relations (discrete string equations) still hold. Therefore, if the double scaling limit exists, 
one should still obtain some solution of the string equation (|j) . 

As already mentioned, it has been shown in ||, within the IDM approach, that this 
is indeed the case, and that there is a one to one correspondence between the weight ratio 
c_|_/c_ and the simply truncated solution of (|j) which is obtained in the double scaling 
limit. To recover this result in the BIPZ approach, let us consider the matrix model integral 

JV 

(P in the large iV limit. The eigenvalue probability density — Xi) becomes the 

i=i 

classical density i>(A). It is convenient to consider the function 

It can be shown, from the saddle point equations for the effective action for v, or through 
the loop equations 0, that F must be of the form 

F(X) = \ [V(X) + VOW] ; Q(X) = V'(X) 2 + 4N(X) (10) 

where iV(A) is some polynomial of degree 

degree N = degreed — 2 = m — 1 (11) 

Generically, Q has 2m complex zeros, which correspond to square root cuts for F. From 
(H) the e.v. density v is proportional to the discontinuity of F along the cuts, and can 
be reconstructed from F. The normalization j dXv = 1 implies that F ~ A -1 at oo and 
this fixes the coefficient of the leading term of N. Requiring that F has only one cut (as 
done for instance in |I(|]) implies that Q has m — 1 double zeros and this fixes uniquely N. 



However, in general F may have several cuts. 

Let us label by a the cuts and by x a the fraction of e.v. along each cut (x a > and 
^2 x a = !)• The e.v. density v must minimize the action 

S = J dXv(X)V(X) - J j dXdiJ,v(X)v(fi) ln\X- n\ + ^r a ^r a - J dXv(X)^j (12) 

where r a are Lagrange multipliers. In fact ( |T0| ) is the most general solution when one 
extremizes S w.r.t. variations of the density which do not change the x a 's. As a conse- 
quence, the effective potential T(A) for one e.v. in the background created by the N — 1 
other e.v.'s 

T(A) = V(X) - 2 +dtiv{n) In(X-fi) (13) 
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is constant along each cut a, and equal to r a 

r(A) = T a if A Get (14) 



so that the total action is 



S = \ E / d\v(\)V(\) + x a T c 



(15) 



The remaining constraints which fix N are the following: 

(a) The e.v. fractions x a must be real 

(b) One must minimize Re(S') w.r.t. variations of the x a 's, subjected to the constraints 
that x a > and that = 1. Since from ( |i~2|) -j^- = T a , this implies that the real 
part of the effective potential is the same along all the cuts: 

Re(r o ) = T unless x a = (16) 

In fact the constraints (a) and (b) can be recast in the same form: 

Im(J a ) =0 , X a = I £-F(\) (17) 

Jc a 2m 

where C a is any contour encircling a pair of zeros of Q. Indeed, if C encircles a cut a, 
I = x a , while if C encircles the end points of two different cuts a and f3, it follows from the 
fact that away from the cuts, T'(A) = V'(X) — 2F(X) (obtained by taking the derivative 
of fll3|) ), that X = ^(1^ — r a ) (see fig. |]). Since there are 2m — 2 independent contours, 
( [LTD fixes the m — 1 remaining coefficients of N. Let us note that if Q has a double zero, 
two independent constraints (|T7J) are automatically satisfied, since for any contour C which 
encircles the double zero and at most one single zero, 1 = 0. Therefore we expect in general 
to have 1 solution of (|17D with no double zeroes, 2m — 1 solutions with one double zero, 
etc. . . Some of these different solutions will be excluded because: 

(c) Some x a 's are < 0. 

(d) The contours of integration for the e.v. cannot pass through one of the cuts a. 
Finally, among the remaining solutions, it is the one with minimal Re(S') (real part of the 
action ([P2"D) which is the physical saddle point. 



Let us discuss explicitly the case of the m = 2 critical point. We start from the 
potential V^(A) = — A 3 /3 + gX. In the critical regime we rescale g c — g ~ a 2 t, A — A c ~ ap 
|TTJ| . t is the renormalized cosmological constant and p the loop momentum. The scaling 



parameter (short-distance cut-off) a is defined so that the double scaling limit is obtained 
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by taking N — > oo, <7 st 2 ing = N 2 a 5 = 1. In the planar scaling limit (N = oo, then a — > 0), 
the general solution ( [TUP for F(A) becomes 

F(X) -> w(p) = a 5/2 -Vp 3 -3tp + c (18) 

where c is the only parameter to be determined in the polynomial N(X) which is relevant 
in the scaling limit. It corresponds to the v.e.v of the puncture operator ||. In the weak 
coupling phase t > 0, the saddle point is the standard one cut solution IflOf ||: 



c = 2t 3 / 2 ; w(p) = a 5/2 \ {t 1 ' 2 -p)Vp + 2 W 2 (19) 

The e.v. are located on the real axis along the cut ] — oo, — 2t 1 / 2 ] (see fig. |3| (a)). The 
action for this solution scales as 

S oc a 5 / 2 t 5 / 2 (20) 

There are other unphysical solutions which violate (c) or (d). 

This solution can be analytically continued into the strong coupling unstable phase 
t < 0. The two complex conjugate solutions describe e.v.'s still located along a single arc 
TT}1 , and correspond to the triply truncated solutions of (01) . From (^D|) they have a purely 



imaginary action. However for t < the constraints fll7| ) have another solution, where 
w(p) has now three branch points. For this solution c is given by 

c = c (-t) 3/2 ; c> (21) 

and it corresponds to the branched distribution of the e.v's depicted on fig. |3] (b) . The 
density of e.v. along the negative real axis vanishes at the real branch point po as y/po — p, 
but there are two arcs starting from po toward the two other complex conjugate branch 
points p±. Therefore an equal fraction x+ = X- of e.v. are sitting along these two arcs. 
Moreover the action for this solution is real and negative, and therefore it is generically 
the dominant one. Away from the negative real axis, the branched solution still exists, but 
with asymmetric branches (x+ 7^ x_), provided that one stays in the sector HI (— n/5 < 
Arg(— t) < 7r/5), and still has a lower action than the perturbative one. Moreover, since 
the constraints (|T7| ) are non-analytic, this solution does not depend analytically of t, (in 
other word, the number c in depends on Arg(t)). 

Finally, in the sectors H and IV, another kind of solutions, with two cuts, is dominant 
(see fig. H (c-d)). These solutions are still non-analytic in t. In the sectors I and V, the 



perturbative analytic solution (|T9| ) is the physical one. 

As is clear from fig. ^|, for the general integration prescription ([7p, one can obtain 
these new, non-analytic, solutions in the sectors H, HI and IV. Since they have lower action 
than the analytic solution, they will dominate the large limit, unless c+ or c_ is zero. 
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This new solution for the e.v. distribution is associated to the simply truncated 
solutions of (J|). This can be seen as follows. First, we have seen that these solutions are 
non analytic in t in the sectors where simply truncated solutions have poles. Since double 
poles of u(t) should correspond to simple zeros of the r function, which corresponds to the 
partition function Z for the matrix model, and since in the large iV limit Z is obtained 
from the action S for the saddle point e.v. configuration by Z = exp(N 2 S), using Cauchy 
formula for the derivative of ln(Z) and Stoke's formula we obtain the estimate for the 
number of zeros in a domain D 

# of zeros / dg dln(Z) _ ^ Ff dgdg d d g ^ _^ ^ 



in domain D J dD 2m dg J J D An dg dg 

Thus the new solution describes a partition function with a positive density of zeros p oc 
ddS in the three sectors H, in, IV, which scales as p oc A^ 2 a 5 / 2 |t| 1 / 2 /(0), where (p = Arg(t). 

One can make the identification more precise by relating the constraints ([17]) to the 
asymptotics of the simply truncated solutions of (||). Following |H| (see |H] jn| for more 
recent discussions) we make the change of variable 

u = t x l 2 V ; T = -t 5 / 4 (23) 

5 

in (|]), which becomes 

U" _ 6C/ 2 + 6 = _V 4 U 

T 25 T 2 y ' 

As \T\ — > oo, U is asymptotic to a Weirstrass elliptic p function Uo(T,Eo), solution of 

{U' Y = 4U% - 12U Q + E (25) 
which is doubly periodic (with a lattice of double poles) with periods 

Oi, 2 = / dU (4U^-12Uo + E )~ 1/2 (26) 
Jc li2 

where are two contours encircling pairs of zeros of the r.h.s. of (|2~5D. In a neighborhood 
of some T = Tq, one can treat Eq as a slowly varying variable Eq(T). From (^) , in the 
local periods coordinates T — To = fiiy 1 + ^2y 2 , Eo varies as 

~ -1 J lt3 (Eo) ; J 1>2 (Eo) = jf dU {4U$- \2U + E ) 1/2 (27) 



Solving the flow equations (p7|) one can check that asymptotically, depends only on the 
argument of T, = Arg(T), but not on its modulus, and is solution of the two constraints 

Re[e ie Ji, 2 (£ )] = (28) 
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But these constraints are exactly equivalent to (|T7|) for ([18]), once we identify E — t 3 / 2 c 
and use the fact that = |Arg(t). 

Before discussing the asymptotics for the loop operators, let us show how the choice 
of weight contours c± in (|7|) fixes uniquely the non-perturbative part of the solution of (f|). 
Let us denote by Z^{0) (resp. F^{0) the partition function ([]]) (resp. it's Logarithm) for 
N e.v.'s with the contour coefficients (§). Taking the derivative of F w.r.t. 9 singles out 
one of the e.v.'s 

dF N 



de 



^-C N [ d\ [ n dX.A^fe-^^ (29) 



where Ci is the contour going from joo to joo. One estimates this integral by first inte- 
grating out the N — 1 last e.v.'s (by using the BIPZ method), in the effective potential 
V(X) = V(X) + jj(V{X) - 21n(A - A)) modified by the first e.v. The resulting effective 
potential for the first e.v. is in general complicated, since it takes into account the back- 
reaction of this e.v. on the bulk N — 1 others e.v.'s. However it takes a simple form if Ao 
is close to the end-point A e of the e.v. distribution, since we obtain 

dFM i f 1 



d6 



where T(A) is the effective potential ([13]) . For the m = 2 critical point, in the scaling regime 
T'(p) = —2w(p), with w(p) given by (|T^). At large iV the integral ( |3T)| ) is dominated by 
the instanton configuration of where the e.v. sits at the top of the potential p = t 1//2 . 



The result, including the contribution of fluctuations around the saddle point, is 

dF i ,o 3~ 3/4 e/o _14~/o + 5/4, „ 1 



1/2 



a «'- i ^«- / 'e-4»^(i + (i)) (3D 

where g s = A^ _1 a _5//2 is the string coupling constant. In the double scaling limit g s = 1, 
and (|3l]) gives the non-perturbative ^-dependence of the string susceptibility u = —F". 
Using the fact that the triply truncated solutions u± correspond to 6* = =f|, one thus 
recovers the results of Q for the non-perturbative part of the simply truncated solutions 
of Painleve I. 

Finally, let us return to the construction of the loop operators. It is very easy to 
check that with the generic choice of contours (0), the finite- N loop equations of the 
matrix model are still satisfied by the loop operators. It remains to understand what 
is the continuum limit of these operators. In string perturbation theory, the operator 
which creates a macroscopic loop with momentum p (conjugate to the loop length £), 
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w(p) = die p£ w(£), can be expressed in terms of matrix elements of the operator Q 
given by (0) ^M. For instance the one- loop correlator is 

,00 1 

(w(p)) = / dx(x\——\x) (32) 

and the problem is to define the resolvent G(x;p) = (x\(p — Q)~ x \x) for the generic real 
solutions of (|j), with poles on the negative real axis. G must satisfy the Gelfand-Dikii 
equation 

-2GG" +G' 2 + 4(p + 2u(x))G 2 = 1 (33) 

and generically G has also double poles at the poles of w. In fact there is a unique 
asymptotic prescription for G(x) in the strong coupling regime x — > —00 which is consistent 
with the contour prescription (0) (defined by the ^-parameter), and the specific it. If we 
perform the rescaling (similar to (|23| )) 

X = \x 5 ' A ■ G = x'^H ; P = x'^p (34) 



in the limit \X\ — > 00 ( p3|) becomes 

-2HH" +H' 2 + 4(P + 2U Q (X))H 2 = 1 + 0(— ) (35) 

X 

where Uq(X) is the elliptic function given by (|23|), which is solution of (P3|). There is a 
unique solution of (|35| ) which is doubly periodic with the same periods Qi^(Eq) than Uq. 
It is given explicitly by 

H(X-P) = - P -^W (36) 
V ; y/4P d - 12P + E 



This, together with ( |32] ) and the constraints ([2Sp which fix Eq, gives the same asymptotic 
expression for the one loop correlator w(p) in the non-perturbative phase t < than 
the expression ( |T5| ) that we have obtained previously through the BIPZ approach. This 
achieves the identification of our large iV non-perturbative solution of the matrix model 
with the real solutions of the string equation (|J). These loop operators will satisfy the loop 
equations ]E[]|J15|P6]|, at variance with the operators constructed only in the perturbative 



phase with the prescription of @ H] [14]. Each loop operator will have a single pole in t 
wherever the string susceptibility has a double pole. This is in fact natural, since the string 
susceptibility is the v.e.v. of two microscopic loops (puncture operator). 

Let us summarize and discuss our results for the case of pure gravity. 

The proposal of to take real combination of complex integration contours in the one 
matrix model to obtain real solutions of the Painleve I string equation for pure 2d gravity 
(c = 0) has been formulated here in the framework of the large iV solution of the matrix 
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model a la BIPZ, i.e. in terms of distribution of eigenvalues. We have shown that the strong 
coupling phase, which corresponds to negative values of the renormalized cosmological 
constant t, and in which the string susceptibility has poles, can be described simply in 
term of splitting of the mean-field distribution of the eigenvalues into two branches at the 
end of the e.v. distribution. 

These two different branches can be viewed as two different topological sectors in 
the integral over the e.v.'s. and the non-perturbative ^-parameter which distinguishes 
the different solutions of the string equation is simply the phase difference between these 
two sectors, which has to be specified in the functional integral. Therefore, at a formal 
level, each 9 defines a ^-vacuum of 2d gravity, as in field theories with topological sectors, 
such as 4d non-Abelian gauge theories or some 2d cr-models. The non-perturbative effects 
associated to this 9 parameter can be estimated by simple instanton methods. 

Finally, we have shown that, for each real solution of the string equation vacuum), 
it is possible to construct in a consistent way observables (loop operators), in such a way 
that the loop equations should be satisfied non-perturbatively. 

Of course, many interesting questions are still open. 

It is clear that one can define non-perturbatively the one matrix model for general po- 
tential, and probably reconstruct by adequate choice of contours the real non-perturbative 
solutions of the unstable even m string equations. We shall discuss below the case of the 
double well potential. Similarly, the same recipe can be applied to the multi-matrix mod- 
els, and used to study the general (p, q) string equations (although for the multi matrix 
models there is no simple picture of the large iV limit in term of e.v. distribution). 

The fact that the loop equations are still valid non-perturbatively in the framework 
discussed here is quite appealing. These equations can be derived from various point of 
views: Dyson-Schwinger equations for the matrix models, Virasoro constraints for the 
KdV hierarchy, recursion relations in 2d topological gravity. This is at variance with other 



schemes which have been proposed for defining non-perturbatively 2d gravity [[19[ |20|| 

One important issue has to be properly understood. In the strong coupling phase 
(t < 0) the partition function Z = exp(F) has zeros, and the loop operators have poles. 
This implies that the non-perturbative real solutions of 2d gravity that we have discussed 
here should suffer from non-perturbative violation of positivity, even in the weak coupling 
regime t > 0. It remains to understand what this really means when one formulates these 
solutions in term of string field theories in low dimensional backgrounds, in particular for 
positivity and unitarity. 



The c = 1 matrix model solution studied in |yj does not suffer from the kind of 
instability of the c = model, since it corresponds to free fermions in an inverted harmonic 
potential, with the two wells of the potential filled at the same Fermi level. Consequently, 
although the string perturbation theory for the c = 1 model is not Borel summable, there 



9 



is a well defined summation prescription which allows to reconstruct this non-perturbative 
solution. Does the kind of ideas discussed here allow to construct other non-perturbative 
solutions of the c = 1 model? 

Finally let us briefly discuss the case of the double well potential 

V(X) = ^A 2 + \\± (37) 

For 11 < large enough, the e.v. are distributed along two cuts (symmetric under A <-> —A). 
At the critical point fx c , the two cuts fuse (at A = 0) into one segment. In the double scaling 
limit the string equation for this critical point is the Painleve II equation ||22|| . Recently, 
Brower et al. showed that by relaxing the parity condition rc n (\) = (—l) n n n (—X) on 
the orthonormal polynomial and on the associated solutions of the recurrence equations, 



new symmetry breaking solutions of the model could be obtained [|12| . This can be easily 
understood by considering the three independent paths of integration for the potential 
(p7[). In addition to the real axis C r we can also integrate over the paths C± going from 
— oo to ±zoo. The most general weight factors for these paths which give a real partition 
function are 

c r = 1 — 2x ; c± = x ± id (38) 

If c± 7^ 0, in the strong coupling phase fx > fx c the e.v distribution is no more the one 
cut solution but a cross-shaped distribution with four cuts meeting at the origin. Setting 
x 7^ breaks explicitly the symmetry A <-> — A and should allow to recover the symmetry 
breaking solutions of |E| (which differ from the standard solution by subdominant terms 
of order 1/N 2 ). Setting x = but 9^0 gives solutions which differ non-perturbatively 
from the standard one, and which correspond to solutions of the Painleve H equation with 
(simple) poles on the negative real axis. 

These considerations can be extended easily to the multicut matrix models studied in 

m 
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Figure Captions 

Fig. 1. The three contours for the cubic potential 

Fig. 2. The contours in (|T7| ) corresponding to constraints (a) and (b) 

Fig. 3. Schematic picture of the e.v. distribution (black line) and the Re(r) > domains 
(grey) in the p complex plane for the generic solution of the m = 2 critical point: 
(a) real t > (and sectors I and V); (b) real t < (and sector in); (c) t in sector 
H; (d) t in sector IV. 
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